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Noise Spectrum of a Diode With a Retarding Field 1 

By Jacob J. Freeman 



A general expression for the noise-power spectrum generated by the random emission 
of electrons of arbitrary trajectories within a waveguide is obtained. A numerical solution 
of the potential distribution within a plane diode is effected, and a family of curves from 
which the potential distribution may be obtained is graphed. A graphical comparison be- 
tween several actual potential distributions and Von Laue's solution is given. The above 
results are utilized to derive the equivalent mean-square fluctuation current due to the space 
charge within a diode for two cases of potential distribution, namely, the linear distribution, 
and that distribution which occurs at the neighborhood of the beginning of the retarding 
field. For the case of the linear potential distribution, the equivalent noise temperature of 
the diode conductance is shown to be equal to the cathode temperature. 



I. Introduction 

A considerable increase at high frequencies in 
both the shot noise and conductance of a nega- 
tively biased diode over the theoretical values 
computed on the basis of zero transit angle has 
been reported in the recent literature by several 
investigators. C. N. Smyth [1] 2 observed ex- 
cessive damping and noise from a negatively 
biased diode at 3,300 Mc/s, and suggested it was 
due to currents induced by space charge. A. Van 
der Ziel [2] and A. Versnel [3] reported preliminary 
measurements on the equivalent noise current, 
and on tin 4 conductance of a negatively biased 
diode at wavelengths of 5.8 and 7.25 m, over a 
range of anode voltages. 

In this paper the noise-power spectrum of a 
diode w T ith a retarding field 3 will be investigated 
theoretically, taking into account only the noise 
power generated by those electrons with .insuffi- 
cient energy to reach the plate, in order to be able 
to evaluate the equivalent mean square current 
fluctuations due to space charge alone.- Because 

1 A dissertation submitted to the faculty of the Graduate School of Arts 
and Sciences of the Catholic University of America in partial fulfillment of 
the requirements for the degree of Doctor of Philosophy. 

2 Figures in brackets indicate the literature references at the end of this 
paper. 

3 By a retarding field is meant one in which the slope of the potential curve 
is always negative. 

4 For brevity, the term "space charge" will henceforth be used to designate 
the* aggregate of electrons that have insufficient energy to reach the plate. 



of the mathematical intractability of the potential 
distribution, the analysis is carried out only for 
the case of large saturation current, when the 
slope of the potential curve at the anode is small, 
and the case of small saturation current, when the 
potential distribution is essentially linear. For 
the latter case, the equivalent temperature of the 
conductance is shown to be equal to the cathode 
temperature. 

Coincidentally, the noise power spectrum gen- 
erated by electrons of prescribed trajectory within 
a waveguide of arbitrary cross section is developed, 
partly because the results can be immediately 
applied to the specific problem of a diode across 
a transmission line, and partly because the results 
are necessary for intended investigations of noise 
from transit-time devices within waveguides. 

II. Noise Spectrum 

In this section 5 will be derived an expression 
for the noise power generated by random emission 
within a waveguide, or a transmission line. The 
latter is to be regarded as a special case of the 
former, in which the particular mode propagated 
is the TEM mode, and the word "guide" will be 
used to designate both. Although in this paper 
only the particular case of a retarding-field diode 

* MKS units are employed throughout this paper. 
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connected to a transmission line will be consid- 
ered in detail, the noise spectrum of a diode 
inserted within a cylindrical guide of arbitrary 
cross section will be treated, since the generaliza- 
tion involves little further complication, and it 
seems desirable to have an explicit treatment of the 
problem for purposes of future reference (e. g., 
computation of the noise power generated by 
electronic transit-time devices within waveguides) . 
It is assumed that the cathode and anode surfaces 
of the diode coincide with parts of the guide sur- 
face, so that there are no discontinuities due to 
the electrode structure. In other words, elec- 
trons emitted from a given portion of the guide 
surface describe a trajectory that is prescribed 
for electrons with given initial velocity by a given 
direct-current potential and impinge on another 
portion of the guide surface. It is assumed that 
the given direct-current potential has no minimum, 
so that the results of this section will apply either 
to the case of a temperature-saturated diode, or 
one with a retarding field, in which the anode 
voltage is sufficiently below the cathode, so that 
the slope of the potential curve at the anode is 
negative. If no direct-current potential minimum 
exists, the perturbations of the trajectory of elec- 
trons with a given initial velocity due to small 
fluctuations in the space-charge potential will be 
small, and it is assumed that they can be neg- 
lected. If a potential minimum exists (fig. 1), 
there is a discontinuous change in trajectory for 
electrons having just enough initial energy to 
reach the potential minimum, so that a small 




fluctuation of potential may cause a large fluctua- 
tion of trajectory for this class of electrons. 

Accordingly, then, it may reasonably be assumed 
that each electron with a given initial velocity 
traverses a fixed trajectory independently of all 
other electrons; thus, one may regard the result- 
ing electromagnetic pulses generated by each 
electron as random events. 

The guide is considered to have its axis in 
the ^-direction, and to be terminated at either 
end by a matched load, with the diode in the 
center. The immediate problem is to compute 
the average noise energy absorbed by each load, 
due to electrons having an arbitrary trajectory. 
Only the case where the range in frequency is such 
that the guide supports only the dominant mode 
will be considered. 

Then if E=E exp(jcot) represents the electric 
field of the dominant mode generated by a current 
density J=J exp(jcot), it is shown in appendix 1 
that 

Eo=-^Hr)fj(r')Mr')exv(-ja\z'-z\)dT'. 

(1) 

Here n is the magnetic permeability, and dr' is an 
element of volume occupied by the current density. 
w=2rrj is the angular frequency, a the phase con- 
stant for the dominant mode, and 6 is the corre- 
sponding normal mode [4]. Some properties of 6 
are briefly discussed in appendix 1. r(x,y,z) and 
r'(x f ,y f ,z') are the coordinate vectors of the field 
and source points, respectively. If the current 
density J(r') has negligible cross section, J(r') 
dr'=I(r')ds', where J(r') is the current flowing 
at (x',y',z'), and ds' is the incremental displace- 
ment along the path of current flow, s', corre- 
sponding to dr'. Equation 1 then becomes 



K 



Figure 1. Voltage distribution of diode with minimum. 



(r)=-^ a b(r)jl(nds^Mryxp[ja(^-z)],z>z\ 

(2) 

When the current I(x',y',z') consists of an electron 
of charge e following the trajectory s', 

I=eb(t-t') = £- f " exp[;'co(*-0]rf«, (3) 

AW J -a, 

where t' represents the time at which the electron 
is at (x', y', z'), and 8(t—t') = 0, for t?*t', and 
J*5(t—t')dt=l. Accordingly, the field generated 
by an electron is 
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JT(0 — ^JJlfe«P !/<-«— MS 

f U v(t').b exp [j(az'-G)t')]dt' y (4) 

where the variable of integration is changed from 
s' to t', v(t')=d8'/dt', and t is the transit time of 
the electron. 
Let 



j(»)= P°i>(0-& exp tf(as'— «*')]#'. ( 5 ) 



Then, 



^(0 = -^J_ ro Jexp tf(«rf-as)]«(«)<fe. (6) 

Since H=j/ijlco curl i?, one finds similarly that 

H=—j- exj)[j(G>t—az)]s(<*))da) 

f — jakXb+ curl b l ,^ 

Since it is desired that E(t) be real, a is defined 
as an odd function of w, so that the imaginary 
parts of the integrals in eq. 6 and 7 vanish, since 
they are odd functions of w. 

The instantaneous power flow due to the pas- 
sage of a single electron, is, by Poynting's 
Theorem, 

P(t)= [ EXH-nda, (8) 

J 8' 

where n is a unit normal to the surface S". 
If one utilizes the orthogonality properties of the 
normal modes, as stated in eq. 94 and 95 in 
appendix 1, namely, 

I b r b h da=b jh , 
and 

J , [bjXiVzbnl.kda^O, 

it is found on substituting eq 6 and 7 in eq 8, that 

P( ^ = lfr7pJ ro a eXP U(^-az))8{o)d<a 

exp [j(o)t— QLz)}s{u>)du. (9) 

Since the pulses of energy are random, or statis- 
tically independent, the average power flow, by 
Campbell's Theorem [5], is 



--P(t) = \ (~ P(t)dt, 

J —CO 



where X is the average number of electrons emitted 
per second, and \e=I, the direct current. 

P =lfi II jl a ^ [J^t-az)H.)d. 

s(w) exp [j(ut-az)]dudt. (10) 

By ParsevaPs Theorem [6], 

r.Mt)MTt)dt= f" F 2 V)F 1 (±f)df, (11) 

J — co J -co 

where /i and f 2 are the Fourier transforms of Fi 
and F 2 . 

Substituting eq 11 in eq 10, 



P=Ji 2 



J 00 



— s(co)s( — co)rfw. 



Or, since s(— co)=s*(co) where s* is the complex 
conjugate of s, and since s*s is even in co, 

p=^n-Kaoi 2 <^. (12) 

If one defines un/a=Z, the wave impedance, then 

P=||f"z(«)|8(«)| 2 d«. (13) 

We consider now a plane parallel transmission 
line, with planes a distance a apart, and of width 
d (fig. 2). For this case, b = i/-y/ad, and a=u^ixe, 
where i is a unit vector in the ^-direction, and € is 
the specific inductive capacity of free space. 
Hence, 



jp el 



el i— f °° 

= I^ /X€ J |s(co)| 2 rfco, (14) 



where 



a( „), MW^Mm*-*) # (15) 

Jo -yjad 

If the trajectory is considered to be in the x,y, 
plane, and the linear dimensions of the cathode 
small compared to the wavelength corresponding 
to any frequency in the range considered, then 
one may take z' = and let the origin of coordi- 
nates be at the cathode. Equation 15 becomes 



Noise Spectrum of Retarding Field Diode 



77 



«(« 



h: 



* p( Q.fexp (-jcot')dt' 
^Jad 



(16) 



In order to compute the power spectrum from 
a diode with a retarding field the trajectories must 
be computed, and this requires that the potential 
distribution for the case of a retarding field be 
determined. This is indicated in the next section. 




Figure 2. Schematic of parallel-plane transmission line. 

III. Potential Within a Plane Parallel 
Diode Subject to a Retarding Field 

1. Rigorous Solution 

The potential distribution within a conducting 
diode has been thoroughly investigated by Fry [7] 
and Langmuir [8], among others. However, Fry 
limited his investigations to the case where the 
potential has a minimum. Von Laue [9] has 
exhaustively treated the potential distribution of 
a diode for the case of a retarding field, where it is 
assumed that the plate curient, I P) makes a negli- 
gible contribution to the potential in comparison 
to the contribution made by those electrons that 
never reach the anode. In other words, Von Laue 
assumes that the electron gas between the two 
electrodes is in statistical equilibrium, and that 
the space-charge distribution is governed by the 
familiar relation, 



p = Po exp 



(~e\) 



where p is the space charge per unit volume. 
Since the rigorous potential distribution for the 
case of a retarding field has never been investi- 
gated, it was believed desirable that this be done, 
in order to check the accuracy and range of Von 
Laue's solution. 

The separation, a, between cathode and anode 
(fig. 3) is considered to be sufficiently small with 
respect to the area of the electrodes so that edge 
effects may be neglected, and the potential 
distribution is governed by 



d 2 V 
dx 2 = 



(17) 



where V is the potential, p the charge density, and 

€ = l/367r.l0~ 9 farad per meter. 

The anode is assumed to be at a voltage V\ 
sufficiently negative with respect to the cathode 



CATHODE 



ANODE 




-V, 



Figure 3. Voltage distribution cf diode with retarding field. 

so that the slope of the potential V(x) is always 
negative. 

The following analysis follows Fry's analysis 
of the a-space precisely, and is outlined here for 
completeness. The only difference, of course, is 
that here the constant of integration is not zero, 
as it is when one has a potential minimum. If 
v Q represents the initial velocity of an electron, 
and v is its velocity at position x, 



1 



~Mv 2 +eV=j:Mv { 



(18) 



where e~ — 1.60 10~ 19 coulomb, is the charge of 
the electron, and ilf=9.10~ 31 , is the mass of the 
electron in kilograms. 

At a given position x, those electrons having 
initial velocities such that lj2Mvl<^eV never 
reach x; those having initial velocities such that 
eV<C\\2Mvl<^eV x pass by x twice, going and 



78 



Journal of Research 



returning. Those having initial velocitis such 
that \j2Mviy>eVi are captured by the plate and 
pass the point only once. Accordingly, the charge 
density is given by 



p(x) = 



) = er ^dv +2ef' 7V '^W (19) 

J jrVi V J -y/rVix) V 

where n(v ) is the number of electrons emitted 
per second per square meter with initial velocities 
between v and v -{-dvo, and r=2e/M. Substituting 
eq 19 and 18 in eq 17, 

d 2 V_ eV f 00 n (v ) dv . C VrTT n (v ) dv 1 f 

Integrating eq 20 between the limits % and a, we get 

fdV\ 2 _(dV\ 2 ^ 
\dxj a \dx) 

2MV f ro , 

~~T n (v Q ) m vl — r V x dv - 

_n(v )^lv 2 —rVdv - 
2 _ n{v Q )^vl-rVdrA' (21) 

J JrV -J 



Since the initial velocities are distributed 
according to the Maxwell-Boltzmann distribution, 



r ^ NM 



X H~2F- 



2kTJ' 



(22) 



in which A 7 is the total number of electrons emitted 
per square meter per second, k is Boltzmann's 
constant, A:=1.38 10~ 23 , and T is the absolute 
temperature. 

Io=Ne, (23) 

where 1 is the saturation current density. If one 
changes to reduced variables, 



and 



where 



£=Vr (*-<*), 

"" kT ' 

^2t<M el p 



(kT) 



3/2 7 



(24) 
(25) 

(25a) 



and /„ is the plate current, then eq 21 reduces, 
when n(v ) is replaced by its value, (eq 22), to the 
following : 
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(3H*)>'< ,+ " , *>-( ,+3 &> < 26 > 

Here erf rj is the error function, 

erf x=-^= Pexp(-/ 2 >/f, (27) 

VirJo 

and (drj/d£)l is the value of (drj/d^)' 2 at {=0, that 
is, at the anode, x=a. 
For brevity, let 



and 



/(ry) = (1 + erf V^7)^-(l + 7- V^)* 



Accordingly, eq 20 may be written as 

Since £ is always negative, 

Jo i/M+m 1 



(28) 



(29) 



(30) 



(31) 



The value of m is determined by the boundary 

condition at the cathode 

Al the cat bode,x=0,f=— ayj—£i, ti^eVijkT=tii, 

so that 



t, 



Jo V/(*)+m 2 



(32) 



fixes the value of w 2 . 



Since for large values of 77, /(t?)->2^, the upper 
limit of the integral in eq 32 may be replaced by 
infinity, so that one has, approximately, 



x 



drj 

(/0?)+m 2 )*' 



Since & is proportional to /|, it may be inferred 
from the above equation that a large value of 
m 2 implies a small value of plate current, and that 
a small value of m 2 implies a relatively large value 
of plate current. 

North's paper [10] tabulates /(r?) (his ^(17)), and 
tabulates the solution of eq 14, for the case 
m 2 =0, (his £ a ). Using the computed values of 
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j(rj) from North's paper, £ was determined as a 
function of 77 by numerical integration, for selected 
values of m 2 , and the results are tabulated in 
table 1 for 0^7^20, and are graphed in figure 4 
for 0^77^10. To determine £ for tj>9, one 
observes that for the tabulated values of m 2 , 
f(r})-\-m?~2e r ', for ??>9. Accordingly for t?>9, 
eq 14 yields 



family of £, rj curves plotted in figure 4, one deter- 
mines the particular curve of the family that 
passes through the point &, 771 . From this curve 
*?(£)>(*? <Ci7i>£<Cfi)> the potential V (x) within 
the diode may be determined from 

V(x)_ 



f(u)=f(9) + V2(e- 4 - s -«-' A ), 



(33) 
larger 



x_ 
a~ 



— y 



Zl 



(34) 
(35) 



from which £(77) may be calculated for 
values of 77 than are tabulated in table 1. 

To determine the potential distribution of a 
given tube, one first determines fr and 771 from the 
given values of cathode temperature, plate cur- 
rent, or saturation current, plate voltage, and 
plate-cathode separation, a. Then, from the 



Table 1. Values of £= — I 



dx 



V/ (x)+m 2 



according to the definitions of £ and 77 (eq 24 and 

25). 

2. Approximate Analytic Expressions for the 
Potential 

In order to obtain an analytic approximation 
for the potential, one must examine approxima- 
tions to j (77) that will make the equation 

as a function of rj, for selected values of m 2 



n 


m2 = 0-£ 


W 2 = .05-£ 


w = 0.1-£ 


w2=0.2-£ 


m2 = 0.3-£ 


W 2=0.5-£ 


W2=l_£ 


ra2=3-£ 


7n2=10-£ 


ra2=20-£ 



































0.01 




0. 0426 


0. 0309 


0. 0221 


0. 0181 


0.0141 










.02 




.0814 


.0601 


.0437 


.0359 


.0280 


0. 0199 


0.0115 






.03 




.1171 


.0880 


.0652 


.0533 


.0418 










.04 




. 1503 


. 1147 


.0858 


.0704 


. 0554 


.0396 


.0230 






.05 


0. 4281 


.1814 


.1402 


. 1057 


.0873 


.0688 










.06 




.2108 


.1648 


.1252 


.1038 


.0821 


.0590 


.0344 






.07 




.2386 


.1884 


.1443 


.1200 


. 0952 










.08 




.2650 


.2111 


.1628 


.1360 


.1082 


.0782 


.0458 






.09 




.2902 


.2337 


.1808 


.1516 


.1210 










.10 


.5941 


.3141 


.2542 


.1984 


16.70 


.1337 


.0971 


.0572 






.20 


.8170 


. 5171 


.434 


.3565 


. 3080 


.2528 


.1881 


.1130 






. 30 


.9785 


.6677 


.578 


.4887 


.4285 


.3590 


.2730 


.1673 






.40 


1. 1081 


.7918 


.697 


.6006 


.5339 


.4545 


.3521 


.2200 






. 60 


1.2173 


.9077 


.801 


.6985 


. 6277 


.5409 


.4258 


.2710 


0. 1545 


0. 1105 


.60 


1.3120 


.9903 


.891 


.7866 


.7116 


. 6195 


.4945 


.3203 






.70 


1. 3956 


1.072 


.972 


. 8638 


.7876 


.6914 


.5586 


.3678 






.80 


1.4704 


1.146 


1.045 


.9348 


.8568 


.7574 


.6184 


.4133 






.90 


1. 5380 


1.213 


1.111 


. 9994 


.9201 


.8184 


.6744 


. 4574 






1.0 


1. 5996 


1.273 


1.171 


1.059 


.9786 


.8756 


.7269 


.4995 


.2997 


.2174 


2.0 


2.0134 


1.714 


1.610 


1.488 


1.402 


1.302 


1.120 


. 8355 


. 5461 


.4106 


3.0 


2. 2338 


1.933 


1.829 


1.701 


1.615 


1.520 


1.331 


1.040 


.7259 


.5654 


4.0 


2. 3615 


2.059 


1.954 


1.837 


1. 741 


1.647 


1.456 


1.164 


.8454 


.6762 


5.0 


2. 4376 


2.134 


2.029 


1. 903 


1.817 


1.722 


1. 531 


1.239 


.9206 


.7490 


6.0 


2. 4834 


2.179 


2. 075 


1.947 


1.864 


1.769 


1.577 


1.282 


.9669 


.7946 


7.0 


2.5112 


2.206 


2.103 


1.976 


1.891 


1. 795 


1.605 


1.312 


.9951 


.8227 


8 


2. 5280 


2.224 


2.119 


1.992 


1.907 


1.812 


1.621 


1.329 


1. 0122 


.8398 


9 


2. 5382 


2.234 


2.129 


2.002 


1.918 


1.822 


1.631 


1.339 


1. 0227 


.8502 


10 


2. 5444 


2. 2402 


2. 1352 


2. 0082 


1. 9242 


1. 8282 


1. 6372 


1. 3452 


1. 0289 


.8564 


12 


2. 5504 


2. 2462 


2.1412 . 


2. 0142 


1. 9302 


1. 8342 


1. 6432 


1.3512 


1. 0349 


.8624 


14 


2. 5526 


2. 2484 


2. 1434 


2. 0164 


1. 9324 


1. 8364 


1. 6454 


1. 3534 


1. 0371 


.8646 


16 


2. 5534 


2. 2492 


2. 1361 


2. 0172 


1.9332 


1. 8372 


1. 6462 


1. 3542 


1. 0379 


.8654 


18 


2. 5537 


2. 2495 


2. 1445 


2. 0175 


1. 9335 


1.8375 


1. 6465 


1. 3545 


1. 0382 


.8657 


20 


2. 5338 


2. 2496 


2. 1446 


2.0176 


1. 9336 


1.8376 


1. 6466 


1. 3546 


1. 0383 


.8658 
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Figure 4. Family of £, rj curves. 



■i 



-r 



cte 



->Jf(x)+m 2 



integrable in terms of tabulated functions. For 
7] large, /(ij)~2c, and as rj-^0, f(rj)^7]. 

For T= 1,000° K r ~=ll.7, and i 7 = 11.7(F 1 -V r ). 

For negative plate voltages of one volt or more, 
r? will be larger than, say 2, over most of the volt- 
age range, so that the approximation f(n)~2 
(e 71 — 1) will be fairly good over most of the range 
of rj. This is illustrated in table 2, where /(r?) 
and 2(e'— 1) are tabulated for 0<7j<9. For 

Table 2. ^4 comparison of f(v) and 2{e r > — l) 



V 


f(v) 


2(6-7-1) 











0.2 


0.294 


0.44 


.4 


.716 


.98 


.6 


1.27 


1.64 


.8 


1.98 


2.46 


1 


2.88 


3.44 


2 


11.8 


12.78 


3 


36.9 


38.18 


4 


105.6 


107 


5 


293 


295 


6 


803 


804 


7 


2,190 


2,191 


8 


5,960 


5,960 


9 


16,200 


16, 204 



large values of rj the approximation is quite good, 
but becomes poor for 0<r?<2. On the other 
hand, the error in £ due to the approximation is 
reduced by the addition of m 2 , especially for 
larger values of m 2 . Accordingly, let 



-h: 



dx 



V2^+m 2 -2 V2J0 -y/e x +A 2 ' 



LP 

/2J0 



dx 



where 



*-?-!• 



Equation 36 is a solution of 
d 2 v_ 



d?~ e "' 



(36) 



(37) 



(38) 



which is Von Laue's expression [9], inferred from 
the assumption that equilibrium exists. It is 
equivalent to 



dW_ /o V 



d.r 



2tM 

kT -It 

- e > 



(39) 



which in turn implies that the space charge 
distribution is 



_eV 
kT 

p=p e , 



where 



Po 



='«V : 



27T.U 

kT' 



(40) 
(41) 



If one integrates the rigorous expression for p in 
eq 19, one gets 



- V °Y kT e 2 



(42) 



Since 



erf 



^jkT 



( Vl -V)=eri Vn«l, 



for r?^>2, it is seen that the approximate expres- 
sion for p approaches the rigorous expression 
fairly closely for r?>2, and accordingly, the 
approximate expression for the potential must 
also approach the true potential in this region. 

For -1^,4 2 <0, (0^m 2 <2), the integral of 
eq 36 is given by 
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-£=- 



V2 



""2 



[cos- 1 (-^l-y^ 2 ) 



COS 



Ml 



For ,4 2 =0, (m 2 -2), 
For.4 2 >0, (m 2 >2), 



(43) 



(44) 



-{= 



"2 

sinn" 1 (V< /2 a/^-I ) 1 (45) 

The value of m 2 in the above equations is fixed 
by the boundary conditions at the cathode, 
namely, v=Vi when £=£i. For large values of 
m, eq 45 gives 






(46) 



and the resulting potential distribution is, of 
course, the same linear function of distance one 
gets in the absence of space charge. 

In order to compare eq 43, 44, and 45, which are 
Von Lane's solutions, with the rigorous solutions 
of eq 32 obtained numerically, values of the two 
solutions are tabulated in table 3, for 7?i=10, and 
for m 2 =0, and ra 2 =3. For m 2 =3, the agreement 
between the two solutions is excellent, and even 
for m 2 =0, the worst possible case for a given value 
of rfi, the agreement is fair. 

Unfortunately, the expressions for the potential 
obtained from eq 43 and 45 do not permit the time 
of transit of an electron as a function of its position 
to be calculated explicitly, which is necessary for 
an evaluation of the noise power spectrum. 

However, the potential distribution given by eq 
44 approximates fairly well the shape of the po- 
tential curve for small m 2 , that is, 0<m 2 <20 
(large / ). In figure 5 graphs are plotted of the 
true reduced potential distribution, V/Vi = 
(vi — v)lvi as a function of the reduced distance 
xja=(£i — £)/£i, for m 2 =0, m 2 =10, and m 2 = 20, for 
?7i=10. Also, the relative potential distribution 
obtained from eq 44 is plotted and is seen to fall 
about midway between the curves m 2 = and 
?w 2 = 20. Accordingly, for the case of large plate 
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Figure 5. Comparison between approximate and 7'igorous 

potential distributions. 

A, m 2 =0; B, approximate formula; C, w 2 =-10; D, w 2 =20; E, m—— W&- 

currents, one may use eq 44 for obtaining an ap- 
proximate potential distribution and thence cal- 
culating the noise spectrum. 

Table 3. A comparison of the reduced potential, 
V\V\ = \ — 77/77!, as a function of the reduced distance, 
x/a=l — £/£i, obtained from the rigorous solution {eqlJ^, 
and from Von Laue's apprcximation (eq 43, 45) 



W2 = 


7W2 = 


= 3 




Rigorous 


Approxi- 
mate 


Rigorous 


Approxi- 
mate 


l-V"i 


l-f/*i 


l-*/*i 


l-^i 


l-^i 

















0. 1 


0. 0025 


0. 0015 


0. 0045 


0. 0048 


2 


. 0064 


.0061 


. 0120 


.0125 


.3 


.0130 


.0138 


.0247 


.0254 


.4 


.0240 


.0264 


.0470 


. 0466 


.5 


.0420 


.0470 


.0789 


.0814 


.6 


0719 


. 0814 


.135 


.139 


. 7 


.1221 


. 1385 


.227 


.233 


.8 


. 209 


.234 


.379 


.386 


.9 


.371 


.412 


.629 


.630 


.95 


.522 


. 566 


.799 


.798 


.98 


. 679 


.719 


.916 


.915 


1.0 


1.000 


1.000 


1.000 


1.000 



For the case of small plate currents, and large 
negative anode voltages, the reduced potential 
distribution approaches a straight line of unit 
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slope, also graphed in figure 5 under the label 
m=— Wfi- Accordingly, the following section 
proceeds with the calculation of the noise-power 
spectrum for the two cases, large plate current 
(m 2 small) and small plate current (m 2 large). 

IV. Noise Power Spectrum for a Diode 

1. Case 1: Noise Power Spectrum for Large Plate 
Current (m 2 small) 

If one substitutes the values of £ and rj from eq 34 
and 35 into eq 44, one gets for the potential 
energy of an electron, U, 



where 



U=Ve=2kTln(l + bx), 



j € 2k 7 1 C 2kr 

a a 



(47) 



(48) 



In this paper only the contribution to the noise 
spectrum made by those electrons that never 
reach the plate, but which return to the cathode 
will be considered, since the effect of those elec- 
trons which constitute the plate current has been 
treated by Rack [11]. Accordingly, one considers 
only electrons whose initial kinetic energy, E, is 
less than eV\. 

For an electron emitted at £=0, the time at 
which it has position x is given by 



h: 



dx 



t/2IM't/E-U(x) 



(49) 



If one substitutes for U(x) its value given by eq 47 
and integrates 



t= 



V 



erf 



irMexp(E/2kT 
kT 26 

E 



\ 2k T' 



D ri V 

ln(l + &a;)"l 



E 



(50) 



for 0^^fo/2, where t is the transit time of the 
electron. 



2i^c=m 

«(«)= f= sin 

-yaa Jo 



m 



so that 



IvM^KkT) . [E 
"\If b~ cri ^2kT 



erf V: 



E 

2k T 



(51) 



(52) 



for 0=S<s=V2. 



According to eq 14 and 16, the noise power, 
P E , absorbed by each load in the frequency range 
A/, due to electrons with initial energy E, is 



_eI E I ix, 

! ~~2 V^ 1 



s(co)| 2 



Aoj 
2t' 



where 



-<«>-r 



v4exp(—jut)dt 
^ad 



(53) 
(54) 



and where 1 E is that portion of the saturation 
current composed of electrons of energy E. 

Since t(x) and t'(x) } the first and second times 
at which the electron passes the point x, respec- 
tively, are related by 



and since 



Letting 



t — 1 t, 
v(t -t)=-v(t), 

^ ) = V^iJo °( ( )exp(-J«<)^- 
J^(*)ex P r>(*-*o)]<2(j- 



(55) 



s(o)) = 



(56) 



Si(u)= I v(t)exp(—joot)dt= exp(—jut)dx 
Jo Jo 



where x m is the maximum value of x, 

1 



(57) 



S(w) = 



^ad 



[^(w) — exp (-M,)s,( — cu)]. (58) 



Substituting eq 52 in eq 57 and 58, 



wt erf <J 2Y£~ ln ( 1 + ^ )X )\ 



2 erf 



v 



2kT 



dx. 



(59) 
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To simplify the above somewhat, the substitution 



2/=y2H 7 ~" ln ( 1+bx )> 



(60) 



is made, and one obtains 



S( " ) = " "^53 Jo Ve Sm 



«f erf 3/ 



2 erf 



v 



E_ 

2kTl 



(61) 



Let 



("V23f27~Jo 



V 



2/cT 



2/6" y2 sin 



cot erf y 



2 erf 



4. 



E 

2k T) 



>dy. 



(62) 



Then 



wi-T-KTrMVsn) <63) 



|«(. 



where the value for 6 given by eq 48 has been sub- 
stituted. 

Since the energy distribution of the electrons is 
Maxwellian, 

7 * = F? exp (- E / kT ) dE > ( 64 ) 

and since the resultant noise power due to un- 
corrected noise sources is simply the sum of the 
powers of the individual sources [5], the resultant 
noise power spectrum due to all the electrons that 
return to the cathode is 

where P represents the noise power absorbed by 
each matched load in the interval AJ=Aoo/2t. 



The diode shunted by two matched transmission 
lines is now replaced by its equivalent circuit, 
namely, a diode shunted by impedance Z /2. The 
equivalent mean square fluctuation current, i 2 is 
defined by 

2P=Z /2i 5 , (66) 

where Z =-yJn/ea/d is the characteristic imped- 
ance of the line. If the impedance of the diode is 
large compared to Z , which is assumed to be true, 
% 2 represents the equivalent mean square fluctua- 
tion current that would flow through a diode short- 
circuited at frequency/. 

From eq 65 and 66, 



326/, 



ffiF<y u f)*j, 



where 



Since 



F{Vuf)= P Vl A\jEl2kT)dE, 



(67) 



(68) 



o)L 



= aco 



2 erf -4E\2kT 
from eq 48 and 51, letting 



, m {E-eVA 



c^au^irMllkT, 

A(^/E/2kT)= f VjE/2/cT \ e -u2 sin j c exp ( ^Z^l \ erf u \ du 

Integrating the above equation by parts, 

A( VE^T)= 2 -^fc^ [f™ cos j c exp [^] erf u ) du- 



(69) 

(70) 



V#/2&Tcos jc exp f^T" 1 ] erf V2?2EP{~|- 
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One observes that A(^E/2kT) = 0, for c— and 
c= oo . Since the expression within brackets is 
always less than 2^Ej2kT in absolute value, 






(72) 



(a) Low-Frequency Approximation 

For a cathode temperature of about 1,()()()°K, 
and with a, the cathode-anode spacing expressed 
in millimeters as a m , one has 



c<5X10- 8 /a ra . 



(74) 



and 



F{VJ) =£' e A^E/2kT)dE< 



2c 2 exp(-eV 1 /kT)' 
(73) 



Since E^eVi, one may approximate the sine of the 
argument in eq 71 by the argument for frequen- 
cies such that c<0.3. 

~E~ 



A(JE/2kT) = c ex P [^j£] £ 

Or, integrating by parts, 



which gives an upper limit to F{V\J) for very high 
frequencies. 



ue~ u2 eriudu. 
(75) 



(76) 



Since A 2 (-y/W/2kT) approaches zero as E/kT and 
for E/kT small, and increases exponentially as 
exrj(E/kT) for E/kT^l, and since one is interested 

in P Vl A 2 dE for an upper limit E/kT=eV l /kT>10, 

the contribution to the integral from large values 
of the argument far overshadows the contribution 
from small values of the argument. Hence, we 
may replace the value of A(^E/2kT) by its 
asypmtotic value for large values of the argument , 



T 2 2 *M , 

r=el P aW pp Aj. 



%' 



(78) 



(The above formula holds only for 
c=aw V^M/4F?<0.3). 
2. Case 2: m 1 Large (/„ Small) 



A(jE/2kT)-- 



c /E-eV, 



2V2 



expl 



"-■eVi \ 
2k T )' 



(77) 



Physically, this means that despite their pre- 
ponderance in numbers, the low-energy electrons, 
because of their small excursions toward the plate 
and consequent short transit times, contribute 
practically nothing toward the low-frequency 
noise spectrum. 
Accordingly, 



x 



eVi £7>2 

A 2 dE=*^-, 



For this case, the potential energy of the elec- 
tron is 

U(x)= j 

v } a 

and the time, /, an electron reaches a point x, is 
given by 

«-VB!lS{f[i-Vr^}««W* 

Following the same procedure as in case 1, one 
finds for the noise power absorbed by each load 
in the frequency interval A/, 



P= 



2k T 



where 



{4/eVl)2 fo 



* Fi ~kT (cos 
Ee \ c^E 



cnjE sin a^E) 2 ,^ , 



(79) 



c=o)^2Ma i /eV 1 . 
The equivalent mean square fluctuation current through the diode is 
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~kT S COS 



vy/E sin <rVffi * , 
-^ -^^dEAj. 



(80) 
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For low enough frequencies, such that 
a-\/eVi<^0.3, the expression in brackets approxi- 
mately equals <x 2 E/9. Since also, 



T- 2 __128Ma 2 oi 2 eI 



\TTJ' 



kT 



4f- 



(82) 



f*eVx 

Jo 



Eh^dE~§(kT)\ for eVJkTyiO, (81) 



For large values of V x , or for high frequencies, 
the upper limit in eq 80 may be taken as infinity 
with negligible error. Accordingly, eq 80 may be 
rewritten as 



kT\eVJ Jo • \ C ^E <?E ) 

In appendix 2 it is shown that 



(83) 



where 



w=a 2 kT, 



(85) 



and F is Rummer's function [12], also known as 
the confluent hypergeometric function. 

Accordingly, the expression for % 2 may be 
written, 



16/qI 



{l + 2w)e- w F(~ 



(I'i^)" 



1 A/. (86) 



V. Equivalent Temperature of the Con- 
ductance Due to Space-Charge 

The shunt conductance, g, of the space-charge 
within a diode has recently been calculated by 
Begovich [13], who finds 6 



4/ e 



kTC 



S^gj (l + 2w)e^F^ |w)-l J, (87) 



try 



1 3 



where w is defined above (eq 85) and F ( o> iv w) is 

the confluent hypergeometric function [12]. 

If the above value of g, and the value of i 2 
given by eq 86 are substituted in the Nyquist 
formula [24], 

T 2 =4kT eQ gAf, 

where T eq is the equivalent temperature of the 
conductance, it is immediately found that 

6 The writer is indebted to N. Begovieh for a personal communication 
amplifying the material contained in reference [13]. 



T, 



(84) 



(88) 



in agreement with the experimental results ob- 
tained by Van der Ziel [3]. 

VI. Summary 

1. A general expression for the noise power 
spectrum generated by the random emission of 
electrons of arbitrary trajectories within a wave- 
guide is developed in section I. 

2. A numerical solution of the potential dis- 
tribution within a plane diode is effected in sec- 
tion II, and a family of curves from which the 
potential distribution may be obtained is graphed 
in figure 4. A graphical comparison between 
several actual potential distributions and their 
approximate representation is given in figure 5. 

3. The results of sections I and II are utilized 
to derive the equivalent mean square fluctuation 
current due to the space charge within a diode 
for two types of a retarding potential distribu- 
tion, namely, the linear distribution, and those 
distributions whose slope at the plate is small. 

4. Utilizing the results of section III, and a 
recently computed value of the conductance due 
to the space charge within a diode, the equivalent 
noise temperature is shown to be equal to the 
cathode temperature. 



The author is greatly indebted to K. F. Herzfeld 
for his many helpful discussions and generous 
guidance during the course of this work. Also, 
the writer is extremely grateful to H. Lyons, chief 
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tion and valuable criticisms. 

VIII. References 

[1] C. N. Smyth, Nature 157, 841 (1946). 

[2] A. Van der Ziel, Nature 159, 675 (1947). 

[3] A. Versnel and A. Van der Ziel, Nature 159, 640 

(1947). 
[4] For properties and applications of normal modes to 
waveguides, see H. Bethe, Formal theory of wave- 
guides of arbitrary cross section, Radiation Labora- 
tory Report 43-26 (Mar. 16, 1943). 
[5] N. R. Campbell and V. J. Francis, J. Inst. Elec. Eng. 

93, 45, (1946). 
[6] G. Campbell and R. Foster, Fourier integrals for 
practical applications, Bell Telephone System 
Monograph B-584, p. 39 (1931). 
[7] T. C. Fry, Phys. Rev. 17, 441 (1921). 
[8] I. Langmuir, Phys. Rev. 21, 419 (1923). 
[9] M. Von Laue, Jahrbuch der Radioaktivitat .mid 
Electronik 15, 204 (1918). 
[10] B. J. Thompson, D. O. North, W. A. Harris, RCA 

Rev. 4, 447 (Apr., 1940). 
[11] A. J. Rack, Bell System Tech. J. 17, 592 (Oct, 1938). 
[12] TV M. MacRobert, Functions of a complex variable", 

p. 346, 3ded. (1947). 
[13] N. Begovich, High frequency loading, Bull. Am. 

Phys. Soc. 23, No. 5 (June 24, 1948). 
[14] H. Nyquist, Thermal agitation of electric charge in 

conductors, Phys. Rev. 32, 110 (1928). 
[15] J. W. Miles, The equivalent circuit for a plane dis- 
continuity in a cylindrical waveguide, Inst. Radio 
Eng. 34, 728 (Oct. 1946). 
[16] T. M. MacRobert, Functions of a complex variable, 
3d ed., p. 145 (Macmillian & Co., London, 1947). 

VIII. Appendix 1 

Let x,y be generalized coordinates suitable for 
describing the guide surface S, and z the usual 
cartesian coordinate which is parallel to the axis 
of the guide. The normal modes, bj of electric 
type for a cylindrical guide are a doubly infinite 
set obtained by two independent processes [15,4]. 

bj=vUjXk, (89) 

where A: is a unit vector in the ^-direction, and 
where 

V 2 ^fe/) + /W^)=0, (90) 

and n-V£//=0 on S, the surface of the guide, n 
is normal to S. 



where \7 2 v j J t-pfD J =(), and Uj=0 on £, except for the 
TEM case for which fo=0, and 



bj=ffii)j(x,y)kTjoLVvj(x,y), 



(9D 7 



nXVfli=0 on S. 



(92) 



(The Pj are, in general, different for tin 1 I \ and 

If the U; and v j are normalized, that is, 

JhU*da=Mda = \, (93) 

where the integral is taken over the transverse 
plane S' of the guide, then the following orthogonal 
properties may be demonstrated [15,4]: 

fs>b r b h da=8 jh , (94) 

where 8 Jh = 0, for jy^h, and §;•*,= 1, for j=h. 

fsibjXcurl b' h )-kda=0 for all j and h. (95) 

The phase constant a corresponding to each bj 
is given by 

a 2 = /xea> 2 -0;, (96) 

and Lj— bj exp ( ± jaz) satisfies the vector wave 
equation 

curl curl Lj=/jl€q) 2 Lj. (97) 

Utilizing the above properties of the normal 
modes, the field within a guide excited by an 
arbitrary current distribution will be determined 
by a method due to Bethe [4]. 

Using Gauss' Theorem, if A and D are two 
arbitrary vectors, 






[4 X curl 2>—2>X curl A] • nda = 

(D.VXVXi-^VXVXfl)(/r, (98) 



where n is the outward drawn normal to S, en- 
closing the volume V whose element is dr. If one 
lets A = E, the electric field which satisfies 



curl curl E—iJL€o: 2 E=—jo)fjiJ y 



(99) 



7 The Pi should not be confused with the letter v used in section I to represent 
velocity. 



where J is the current density, and lets D be a 
solution of 

VXVXZ>-/x€co 2 Z)=0, 

then eq 98 becomes 

\ (EXVXD— DXVXE)-nda = — njw J J • Ddr. 

(100) 



Noise Spectrum of Retarding Field Diode 



87 



Consider a current element within the guide at 
0=2', and take matched loads at the two planes, 
2=0 and z=Zi. Expanding the field E in terms 
of the normal modes, let 



E=y^Aibi exp (—jajz)z^>z f 
3 

E=^2Cjbj exp (+jajZ)z<C.z' 



(101) 



Since E is continuous at z=z'j and since the bj 
are orthogonal, 

A j exp (-jajz') = Cj exp (ja,z') . (102) 

If one substitutes the value of E given by eq 101 
in eq 100, and lets D=bj exp(—jajz), and inte- 
grates the left side of eq 100 over the plane z=0, 
one gets 

- f T^CnlbnXVXibjexpi-jajz)) 

Jz=0 h 



-bjXVX(b h exp(+ja h z)] • kda, 



or, 



-JZCnf[bnXWXb j -ja j b h X{kXb j 

h 

- bjX V X b h -ja h bjX (kX b h )] • kda. 

Because of the orthogonality relation (eq 96), 
the above simplifies to 

— S C* [—jajfb h 'bjda—ja h fb j -b h da=2C j a j j. 

h 

The surface integral over the plane z=z x vanishes, 
and one has, from eq 100, 

^J=-^ f J bjexp(-jajz')dT. (103) 

Substituting eq 103 and eq 102 in eq 101, one 
finally gets 



E= ^^jJ-hexpi-faM'-zl ) dr 

2 j a* 



(104) 



IX. Appendix 2 



Let 



i= (he" kT {cos W^Lsin fVffi'dg, (104) 

Jo ( a-yJE <rE ) 

*) 2 _ 2 £ ( ._ lV -(2^- 2 (iz^), (105) 



cos x sin 



which may be verified by performing the square, 
expressing the result in trigonometric functions of 
2x, and then substituting for each trigonometric 
function its corresponding infinite series and 
combining. 

If we substitute, eq 105 in eq 104, we obtain 

/=(*2Tg(-i)'2 w -v*zy- 1 ^ 1 ! ) o-+i) 2) ' 

(106) 

where we have utilized the formula for the gamma 
function, 



V- 1 



e ax dx= 



(107) 



Using the duplication formula [16] for the y- 
function, 

. (2j)!=2i 7 (2j)-^70 , + l)7(i+^) 



eq 106 simplifies to 



r-*F> 






(108) 



From MacRobert's definition of the confluent 
hypergeometric function, it follows that 



F( « >P;W) =%) Sr j^(«+») 



7(«) n=oy(n+l)y(p+n) 
Accordingly, eq A 109 may be written 

(kT) s < 



(109) 



-2p(i| -«,) +j f(i|,- W )J, (110) 



where w=a 2 kT 
However, since 
F(l,p-l,-w) = l-^=^wF(l,p-w), (111) 



7(p) 



(wy 



j(l + 2 W )^(l|- W )-lJ. (112) 



Finally since 

F{a,p,W) = e w F(p — a.p,— w), 
r(kTf( 



(l+2 W )e-»2f(||^)J. 



(113) 
(114) 



Washington, October 6, 1948. 
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